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The color number A'^c-dependence of the interplay between quark-antiquark con- 
densates {qq) and diquark condensates {qq) in vacuum in two-flavor four-fermion 
interaction models is researched. The results show that the Gs-Hs (the coupling 
constant of scalar (gg)^-scalar (qq)'^ channel) phase diagrams will be qualitatively 
consistent with the case of A'^c = 3 as Nc varies in 4D Nambu-Jona Lasinio model 
and 2D Gross-Neveu (GN) model. However, in 3D GN model, the behavior of the 
Gs-Hp (the coupling constant of pseudoscalar {qq)'^ channel) phase diagram will 
obviously depend on Nc- The known characteristic that a 3D GN model has not the 
coexistence phase of the condensates (qq) and (qq) is proven to appear only in the 
case of Nc < 4. In all the models, the regions occupied by the phases containing the 
diquark condensates (qq) in corresponding phase diagrams will gradually decrease 
as Nc grows up and finally go to zero if A'c ^ oo, i.e. in this limit only the pure (qq) 
phase could exist. 
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I. INTRODUCTION 

The four-fermion interaction models including 4D Nambu-Jona-Lasinio (NJL) model 
l|, 3D and 2D Gross-Neveu (GN) model 0] are not only ideal theoretical laboratories 



to research dynamical symmetry breaking [l|, and in recent years, the color 
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superconductivity 0, S] , but also can find many practical usefulness to Particle Physics 
and Condensed Matter Physics. Hence it is certainly interesting to conduct some more 
deepgoing exploration for these field theory models themselves. 

Let q denote the fermion (called as quark) field and q antiquark field. Since in any four- 
fermion interaction model, the^ couplings of (gg)^-form and (gg)^-form can always coexist 
via the Fierz transformation [8], so once a sufficiently strong coupling of (gg)^-form leads 
to formation of the diquark condensates (gg), then there will certainly exist interplay 
between the quark-antiquark condensates {qq) and the diquark condensates {qq) in the 



ground state of the system, even in vacuum [9|, [lO|, lUl . In vacuum condition, despite of 



absence of net quarks, however, based on the relativistic quantum field theory, it is pos- 
sible that the diquark condensates {qq) and di-antiquark condensates {qq) are generated 
simultaneously. 

Motivated on the above idea, we have researched interplay between the condensates 
{qq) and {qq) in vacuum in two-flavor 4D NJL Q, 2D GN and 3D GN 3 models and 
find that the ratio between the coupling constants of the (gg)^-form and the (gg)^-form 
plays a key role in determining the property of the ground state. It should be indicated 
that some properties of the system in the vacuum could be maintained to the case of finite 
temperature and density. For instance, it has been shown that pj], in a 3D GN model 
with the quark's color number N^. = 3, the conclusion derived from the vacuum condition 
that the model does not contain the phase with coexistence of the condensates {qq) and 
{qq) will be kept to the case of finite temperature and finite quark chemical potential. 
Therefore, the property of a model in vacuum could reflect its some essential characteristic 
and this fact will make the research on interplay between the two condensates in vacuum 
be given more theoretical significance. In addition, a technical advantage is that all the 
discussions can be made in an analytical way in this case. 

In our preceding work, the color number is fixed to be 3, consistent with phe- 
nomenology of Quantum Chromodynamics (QCD). In this paper, we will extend the 
research to more general case of any Nc- In fact, the large Nc behavior of a four-fermion 
interaction model is always an interesting problem to deserve to be followed. 
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3,0] 



and 



We will still use the effective potential approach employed in Refs. 
obtain the expressions of relevant effective potentials by similar procedures but omit the 
details of these procedures. The mean field approximation will be taken. In this case we 
will use the Fierz transformed four-fermion couplings in the Hartree approximation so as 
to avoid double counting Q|. In selecting the couplings of the (gg)^-form, we will always 
simulate the SUc{Nc) gauge interactions, where two fermions are attractive each other in 
the antisymmetric Nc{Nc — l)/2-plet. The main results will also be shown in some figures 
more clearly. 
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II. 4D NAMBU-JONA-LASINIO MODEL WITH ANY 

With 2 flavor and color massless fermions, the Lagrangian 

£ = qiYdf^q + Gsliqqf + {qil5rq)^] 

+Hs ^{qil->T2\Aq^){fil^r2\Aq), (1) 

where the fermion fields q are in the doublet of SUf{2) and the A^^^-plet of SUc{Nc), i.e. 

g= I ""M ^ = l,■■■,N,, (2) 




g"-^ is the charge conjugate of q and r = (ri, r2, r^,) are the Pauli matrices acting in two- 
flavor space. The matrices Xa run over all the antisymmetric generators of SUc{Nc), Gs 
and Hs are respectively the coupling constants in scalar (gg)^ channel and scalar color 
Nc{Nc - l)/2-plet (gg)2 channel. 

Assume that the four-fermion interactions can lead to the scalar condensates 

{qq) = (3) 

with all the A^^^ color fermion entering them, and the scalar color Nc{Nc — l)/2-plet 
difermion and di-antifermion condensates (after a global SUc{Nc) transformation) 

{q'^il5r2X2q) = S, {qi'^->T2X2q^) = S*, (4) 

with only two color fermions enter them. The corresponding symmetry breaking is that 
SUfL{2)®SUfR{2) SUf{2), SUc{N,) SUc{2), and a "rotated" electric charge Uq{1) 
and a "rotated" quark number Ug{l) leave unbroken [8|. It should be indicated that in 
the case of vacuum, the Goldstone bosons induced by spontaneous breaking of SUc{Nc) 



could be some combinations of difermions and di-antifermions 121]. 
Deflne that 

a = -2G50, A = -2Hs5, A* = -2Hs5*. (5) 
With a 4D Euclidean momentum cutoff A, we can obtain the relativistic effective potential 

H 



y4(^,|A|) = -— + 



AGs 4:Hs 47r2 



{N^a' + 2\A\')A' - (iV, - 2)^ ( In ^ + i 

2 \ a'' 2 



-(a' + \A\'?{\n ^^\,^ + ^ 



(t2 + |A|2 ' 2 



(6) 
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UN, 



y=2x/Nc 




x=HsA^/ Ji^ 
R: y=x/[l+(Ne-2)x] 



Fig.l: Gs-Hs phase diagram of a 4D Nambu-Jona-Lasinio model with the denotations that 
(0, Ai)-pure (qq) phase, (ct2, A2)-mixed phase with (qq) and (qq) and (cri,0)-pure (qq) phase. 

where the conditions ^ and ^ a^ + A^ have been assumed. Through an analytic 
calculations, we have found that the ground states of the system, i.e. the minimum points 
of V^ia, |A|), will be at 



(^,|A|) 



(0, AO 

(fX2, A2) 

(^1, 0) 



if < 



HsA^n'^ > 1/2, 
l/[l + {N,-2)HsAV7r'[ 



< Gs/Hs <1/[1' 

< Gs/Hs < 



(7) 



Eq.(7) gives the phase diagram Fig.l of the 4D NJL model. 

The results qualitatively coincides with the ones obtained in Ref.[l2| except that the 
color number is changed from 3 to N^. Now the pure (gq) phase happens only if the ratio 



Gs/Hs > 2/Nc. This conclusion was obtained in Ref.|5| but no Gs — Hs phase diagram 
was given there. The 2/Nc is just the ratio of the color numbers of the quarks entering 
the diquark condensates and the quark-antiquark condensates. The diquark condensates 
(qq) could emergy only in the region where HsA'^/n'^ > 1/2 and Gs/Hs < 2/Nc. It is 
seen from Fig.l that as Gs/Hs decreases (or Hs increases), we first come to a phase with 
coexistence of the quark-antiquark condensates and the diquark condensates, then to a 
pure diquark condensate phase, the boundary between the two phases is determined by 
the curve Gs = Hs/[l + {Nc — 2)HsA'^ /tt^]. On the other hand, for a fixed Gs and Hs, as 
Nc grows up, the slopes of the lines Gs = 2Hs/N^ and Gs = Hs/[l + (A'c - 2)HsA^/ti'^] 
will decrease. This implies that the regions with the diquark condensates will become 
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smaller and smaller, and finally reduce to zeroes if oo. In the latter limit, only pure 

quark- antiquark condensates could exist. 



1 ^ 




f l\ 




U -1 J 




= -C, 75 
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III. 2D GROSS-NEVEU MODEL WITH ANY A^^ 

The Lagrangian of the model is expressed by 

C = qi^^d^q + Gs[{qqf + {q^l^rqf] 

+Hs{qil-,TsXAq^){.f'i'l5rs\Aq), (8) 

All the denotations are the same as ones in 4D NJL model, except that in 2D space-time 



and Ts = (tq = 1,Ti,t^) are flavor-triplet symmetric matrices. It is indicated that the 
product matrix C'-yi^TsXA is antisymmetric. 

Assume that the four-fermion interactions could lead to the scalar quark- antiquark 
condensates 

{qq) = 0, (9) 
which will break the discrete symmetries 

Xd ■ qit,x) ^ -f5q{t,x), 

Vi : q{t,x) ^ 7^g(t, -x), 

and that the coupling with Hs can lead to the scalar color Nc{Nc — l)/2-plet difermion 
condensates and the scalar color anti- Nc{Nc — l)/2-plet di-antifermion condensates (after 
a global transformation in flavor and color space) 

(g^z75l/A2g) = 6, {qil^lf\2q^) = S* (10) 

which wil 



break discrete symmetries (center of SUc{Nc)) down to Z|, besides xd 
and Vi ISj- Noting that in a 2D model, no breaking of continuous symmetry needs to be 
considered on the basis of Mermin-Wagner-Coleman theorem 16| . 



The model is renormalizable. We may derive the renormalized effective potential (l3| 
containing the condensates (9) and (10). In the space-time dimension regularization 
approach, we can write down the renormalized C in D = 2 — 2e dimension space-time by 
the replacements 

Gs ^ GsM^-'^Zg, Hs ^ HsM^-'^Zh, 
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with the scale parameter M, the renormahzation constants Zq and Z^. In addition, the 
7^^ in L will become 2^/^ x 2-^/^ matrices. 
Define the order parameters 

o = -2GsM^-^Zg(I), a = -2HsM^-^Zh6, (11) 

which will be finite if Zq and Zh are selected so as to cancel the UV divergences in and 
6. In the minimal substraction scheme, 

Zg = 1 Zh = 1 --. (12) 

Ti e TV e 

The corresponding renormalized effective potential in the mean field approximation up to 
one-loop order becomes 



V2{cr,\A\) 



a 
2^ 



2 In 



+ (iV,-2)ln— + iV, 



|A|2 



AH. 



|A| 



7r 



In 



M2 



|A| 



+ 1 



27re~^M2, 



(13) 



where 7 is the Euler constant. V2(o", |A|) contains the two order parameters a and |A|. 
By almost identical demonstration to the one made in Ref. 13| (the only change is that 
the color number of quark varies from 3 to A^^c), we obtain that the ground states of the 
system i.e. the minimal points of V2(cr, |A|) will be at 



(^,|A|) 



(0, AO 
(0^2, A2) 
(^1, 0) 



if < 



Gs/Hs = 
< Gs/Hs < 
Gs/Hs > 2/N, 



(14) 



The results (14) can be drawn in Fig. 2 which is the Gs — Hs phase diagram of the 2D 
GN model. 



The phase diagram are similar to the one obtained in the case of A'^c = 3, including 
that formations of the condensates do not call for the coupling constant Gs and Hs 
having some lower bounds, however, the boundary between the pure quark-antiquark 
condensate phase and the coexistence phase with quark-antiquark condensate and diquark 
condensates is now replaced by Gs/Hs = 2/Nc. It is seen that as Nc increases, the area 
in Fig. 2 occupied by the coexistence phase of quark-antiquark condensate and diquark 
condensates will gradually decrease and finally disappear if Nc ^ 00. It is indicated that 
if Nc — > 00, the region corresponding to the pure diquark condensates in the diagram 
(the transverse axis of Fig. 2) will also not exist. In fact, the pure diquark condensate 
phase corresponds to the minimal point of V2(o", |A|) (a, |A|) = (0, Ai), where Ai is the 
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y=Gs/ ^ 

y=2x/Nc 




x=Hs/Ji (0, Ai) 



Fig. 2: Gs — Hs phase diagram of a 2D Gross-Neveu model with the denotations that (0, Ai) — pure 
{qq) phase, (cr2, A2) — mixed phase with {qq) and {qq) and (cri,0) — pure {qq) phase. 



non-zero solution of the equation 9V2((T, |A|)/c}|A| = with (7 = 0. In this case, it is easy 
to find that the determinant of the second derivatives of V2 



A2 B2 

B2 C2 



with A2 = 9^2/9(72, B2 = 9V2/9(7/9|A| = 9V2/9|A|/a(7 and C2 = aV2/9|Ap. 
the solution (0, Ai), we get that 



For 



N, 



2G, 



2Hs 



.-2, M2 
In^r- 



+ 



2(iVc - 2) 



TV 



K2 = 4A2/7r. 

If Nc is finite, then Gs — ^ could make A2 > and K2 > Q thus the solution (0, Ai) 
becomes a minimal point of V2. However, if A^^ OO) then the third term of A2 will have 
the —00^ behavior, hence even if 6*5 = then it is impossible to have A2 > 0. As a result, 
the solution (0, Ai) is no longer a extreme point of V2. This will imply that the region of 
the pure diquark condensate phase corresponding to the transverse axis of Fig. 2 will not 
exist for finite Hs, Gs — and 00. 



8 



IV. 3D GROSS-NEVEU MODEL WITH ANY 

The Lagrangian of the model is expressed by 

C = qzYd.q + Gsliqqy + iqrqr] 

+HpJ2iQr2\Aq'')ifT2\Aq), (15) 

where 7'^(/i = 0, 1, 2) are taken to be 2 x 2 matrices 

'=(::)■-(:;)'■=(:;)- 

It is noted that the product matrix Ct2Xa is antisymmetric, and since without the 
"75" matrix, the only possible color Nc{Nc — l)/2-plet difermion interaction chan- 
nel is pseudoscalar one. The condensates (qq) will break time reversal symmetry 
T : q{t,x) — > 7^g(— t,x), special parity Vi : g(t,x^,x^) 7^g(t, — x^, x^), special par- 
ity V2 '■ q{t,x^,x^) — > 7^g(t, X"^, — x^). The difermion condensates {(fT2X2q) (after a 
global rotation in the color space) will break SUc{Nc) — > SUc{2) and leave a "rotated" 
electrical charge Uq{1) and a "rotated" fermion number U'^{1) unbroken. It also breaks 
parity V : q{t, x) — > J^qit, —x) and this shows pseudoscalar feature of the difermion 
condensates. 

Define the order parameters in the 3D GN model 

cj = ~2Gs{qq), A = -2Hp{fr2\2q), (16) 

The effective potential in the mean field approximation can be expressed by [3] 

^2 |A|2 1 



a 



AGs AHp TT"^ 



3 



+^ [Qa^\A\ + 2|A|3 + {N, - 2)a^ 

+2e{a - \A\){a - \A\f] , (17) 

where A3 is a 3D Euclidean momentum cutoff and the assumptions A3 ^ W ~ |A||, 
A3 3> 0" + |A| and A3 ^ cr have been made. The ground states of the system correspond 
to the least value points of V3(cr, | A|). Since the phase structure of the 3D GN model will 
have qualitatively different A^^c-dependence from the one of the 4D and 2D models, we will 
make more detailed discussions of this problem. 

The extreme value conditions dVs{cr, \A\)/d(j = and dV3{a, |A|)/9|A| = will re- 
spectively become 

/ 1 _ 2N,A^ , ^ , (iVc - 2)a 

^\2Gs 7T IT 

+ -^(a-|A|)(a-|A|)2 = 0, 

TT 
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lAI 



2Hp 7r2 



4As ^ 2|A| 



TT 



+-[a'-^((7-|A|)((7-|A|)^ 



0. 



TT 



(19) 



Define the determinant from the second derivatives of V3(cr, |A| 



K 



A B 
B C 



AC - S^ 



where 



A 
B 
C 



1 2iVeA3 41 A| 2(Ar^-2)a 4^, 

+ ^ + ^^ ^ + -^((7- |A|)((7- |A|), 



9(7^ 26*5 TT^ TT TT TT 

^^Vq 4 

9Vs 



1 4A3 4|A| 4^, 

+ ^ + -^(cT - |A|)(a - |A|). 

TT TT 



(20) 



9|A|2 2i/p 7r2 

The equations (18) and (19) have the following four different solutions which will be 
discussed in order. 

(i) (cr, |A|) = (0,0). It is a maximum point of ^3(17, | A|), since in this case we have 

A = I/2G5 - 2AreA3A' < 0, 
K ^ A {l/2Hp - 4A3/7r^) > 0, 

assuming Eqs. (18) and (19) have solutions of non-zero a and |A|. 

(ii) (cr, I A|)=(cri,0), where the non-zero ai satisfies the equation 



I/2G5 - 2N,A3/n^ + N,a,/7T = 0, 



(21) 



under the condition that G^As/tt^ > l/ANc- For this solution, by means of Eq.(21) we 
obtain 



A = iVc(^i/7r, 
K ^ A' ^ 



1 



2Hp NcGs TT 



+ ^1>0, if 



Hp Nr 



Hence ((Ti,0) is a minimum point of Vs{a, |A|) when Gg/Hp > 2/Nc 
(in) ((7, |A|)= (0, Ai), where non-zero Ai obeys the equation 



l/2Hp - 4A3/7r2 + 2Ai/7r = 0, 



(22) 
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under the condition that HpAs/n'^ > 1/8. For this solution, by using Eq.(22), we get 

1 2JVA ^ 4A,_ 



2Gs TT 

K = 2AiA/tt. (25) 

(a) When < 4, from Eqs.(23) and (25) we have 

.>o.x>oto.|<A^>i 

(b) When Nc > 4, from Eqs.(24) and (25) we have 

A > 0, X > for — — < — -, — > -. 

7r2 ANc 7r2 8 

Noting that in this case the constraint Gg/Hp < 2/A^c is also implied owing to that 

1/SHp < As/tt^ < l/4NcGs. 

Thus in the cases of both (a) and (b), the solution (u, | A|) = (0, Ai) is a minimum value 
point oiVaia, |A|). 

(iv) (cr, |A|) = ((72, A2). In view of the function 9{a - |A|) in Eqs.(18) and (20), we 
have to consider the two cases of (72 > A2 and (72 < A2 respectively, 
(a) (72 > A2. In this case Eqs.(18) and (19) will lead to 

1 TV, _ N^al - 2Ai ^ ^ Gs 2 



> thus — ^ < 



2Gs 4Hp 7r(72 Hp 

since iV, > 2. Then we obtain from Eq.(20) that 

A = ^ > and K = ^ < 0. 

7rcr2 TT^ 

Hence the solution {(J2, A2) is not an extreme value point of V^la, | A|) if a2 > A2. 
(b) £72 < A2. In this case Eqs.(18) and (19) become 

1 2iV,A3 ^ 4A2 + (AT, - 2)(72 ^ ^26) 



2Gs TT 

and 



2Hp TT ttA 

which may lead to the equations 



1 4A3 2A2 2(7? ^ 

+ — ^ + = (27) 



2i/f 



KG 



c^S 



nN^A2 
2(4 - A^, 
7rAreA2 



[(4 - AT,) A^ + (AT, - 2)a2A2 - Ar,a2'] 



A2 + 



A^c-2 
2(4 -AT,)" 



(X2 + 



4N,{Nc - 4) - (ATe - 2f 
4(4 - Ar,)2 



(7i(28) 



11 



and Eq.(20) will lead to that 



A 




X = A.-Ili±l^^^^fli±i., (29) 

Obviously, for Nc > 2, K > will be the condition in which the solution ((72, A2) becomes 
a minimum value point of ^3(0-, |A|). By means of Eqs.(28) and (29) we may discuss the 
relation between Gs and Hp when the ((J2, A2) becomes a possible minimum value point 
of ^3(0-, |A|). The result is obviously A^c-dependent. 

(bl) Nc — 2. In this case A — 0, K — — 16(7|/7r^ < 0, so ((T2, A2) is not an extreme 
point of V3. 

(b2) Nc = 3. In this case, K > or K > become A2 > (VT7 + 4)(T2 or A| > 
(7| + 8(72 A2. Substituting these inequahties into Eq.(28) we get 

11 2 



2Hp NcGs 37rA2 



-(A^ + (72 A2 - 3(7^) 



3(71 . /tt: . . 2, 



(b3) Nc — 4. In this case, K > becomes A2 > {V^ + 2)(72, then we have 

= ^-(A2 - 2(72) > - — - > 0. 

2Hp N^Gs nA2^ ^ ^' TrAa 

The results in (b2) and (b3) indicate that when = 3,4 and Gg/Hp > 2/Nc, the non- 
zero solution ((72, A2) of Eqs.(26) and (27) could be a minimum point of V3((7, |A|). 

we have obtained in (ii) that when Gg/Hp > 2/Nc, the non-zero solution ((7i,0) of 
Eqs.(18) and (19) is also a minimum point of V^la, |A|) and the result is iVc- independent. 
In this case we must determine which one of the solutions ((72, A2) and ((7i, 0) is the least 
value point of Vsi^cr, |A|). For this purpose, we will compare the value of V3((7, |A|) at the 
point ((7i,0) with the one at the point ((72, A2). In fact, by using Eqs.(17) and (21), we 
may obtain that 

V3{ai,0) ^ -Ncaf/GTT (30) 
and by using Eqs.(17),(26)and (27), we can get that 

1 / — 2 \ 

V^3(a2,A2)U^<^^ = -- f A3 + 3A2(72' + ^^(72M . (31) 
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In addition, comparing Eq.(21) with Eq.(26) with Gs and Hp being fixed, we will be led 
to the relation that 

N.ai = 4A2 + (TVe - 2)^2 (32) 
By means of Eq.(32) wc may calculate that 

. f [i(23Ai - 4ai) + |a2A2(8A2 - 7^2)] , if TV^ = 3 
Vs{a,,0)-Vs{a2,A2) = -—I 

[ [Al - faf + |a2A2(2A2 - ^2)] , if iV^ = 4 

< for CT2 < A2 (33) 

This shows that when A^c = 3,4 and Gs/Hp > 2/Nc, the least value point of Vs{a, |A|) 
is ((Ti,0) rather than ((J2, A2). The same conclusion is also true for the case of A^c = 2 
owing to the fact that ((T2, A2) is now not an extreme value point of V^la, | A|), as shown 
in (bl). These results indicate that, when Gs/Hp > 2/Nc, only the pure quark- antiquark 
condensates could exist in the ground state of the 3D GN model. 

(b4) Nc > 4. In this case, the solution ((72, A2) which satisfies K > will be a minimum 
value point of V3{a, | A|). Noting that the right-handed side of Eq.(28) is always less than 
zero when Nc> 4: and this implies that 

Gs/Hp < 2/7Ve when > 4, 

hence we may affirm that if the solution ((T2, A2) is a minimum value point of V3 (cr, |A|) 
when A^c > 4, then the condition Gs/Hp <2/Nc must be obeyed. 

Now the least value points of y3((7, |A|) can be summed up as follows: 



(^,|A| 




if Gs/Hp < 2/N„ HpA^/n^ > 1/8, 

if Gs/Hp < 2/N„ HpA^/n'' > 1/8, GsAs/tt' 



< l/4A^c, 
> l/4A^c, 



for Nr<4: 



for A^c > 4 



if Gs/Hp > 2/N,, 



GsA^/tt^ > 1 /4N^, for all 
(34) 



Eq.(34) gives the Gs — Hp phase diagrams Fig.3 and Fig.4 of the 3D GN model. 

It can be seen from Eq.(34) and Figs. 3-4 that, different from the 4D NJL and 2D 
GN model, the qualitative feature of fermionic condensates appearing in the ground state 
obviously depends on Nc- When the necessary conditions G5A3 > n'^/ANc and Hp A3 > 
7r^/8 are satisfied, the ground state of the system could be in a pure quark- antiquark 
condensate phase only if Gs/Hp > 2/Nc for all a similar result to the one obtained in 
the 4D and 2D model. However, if Gs/Hp < 2/Nc, then whether the coexistence phase 
with the condendates (gg) and (gg) could appear will depend on Nc- The coexistence 
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1/4N, 





y=2x/Nc 

\ 


(oi, 0) 




^ (0, Ai) 





x=HpA3/ Ji^ 



Fig. 3: Gs-Hp phase diagram of a 3D Gross-Neveu model when Nc< 4: with the denotations that 
(0, Ai)-pure {qq) phase and (cri,0)-pure {qq) phase. 



y=GsA3/3r^ 



1/4N, 




1/8 



(0,AJ 
x=HpA3/ Ji^ 



Fig. 4: Gs-Hp phase diagram of a 3D Gross-Neveu model when Nc> 5 with the denotations that 
(0, Ai)-pure {qq) phase, (0-2, A2)-mixed phase with {qq) and {qq) and (cri,0)-pure {qq) phase. 

phase could not exist if iVc < 4 and could do only if iVc > 5 in a 3D GN model. In 
addition, also similar to the 4D and 2D model, as increases, the regions in the phase 
diagram occupied by the phases with the diquark condensates including the pure diquark 
condensate phase and the coexistence phase with quark-antiquark and diquark condensate 
will decrease gradually, and finally go to zeroes as Nc ^ 00. 



V. CONCLUSIONS 



We have researched interplay between the quark-antiquark and the diquark conden- 
sates in vacuum in two-flavor four-fermion interaction models with any color number Nc- 
The results are quahtatively consistent with the ones in A^c = 3 case in 4D NJL and 
2D GN model, however different from the ones in A^c = 3 case in 3D GN model. It has 
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been found that if the ratio of the four-fermion couphngs of (gg)^-form and (gg)^-form 
Gs/Hs (or Gs/Hp) > 2/Nc, the ratio of the color numbers of quarks entering the con- 
densates (qq) and {qq), (and also with sufficiently large Gs in 4D and 3D model), then 
only the pure (qq) phase may exist and this conclusion is A^^c- independent. Below 2/Nc, 
(and also with sufficiently large Hs or Hp in 4D or 3D model), one will always first have 
a mixed phase with the condensates {qq) and (qq), then a pure (qq) phase, except in the 
3D GN model where the mixed phase does not appear if Nc < 4, the latter includes the 
known result in A^^^ = 3 case. A common characteristic in all the models is that as 
increases, the phases with the diquark condensates (qq) including the pure (qq) phase 
and the mixed phase with the condensates (qq) and (qq) will occupy smaller and smaller 
region in corresponding phase diagram and finally disappear as A"c oo. In that limit 
only the quark-antiquark condensate phase become possible. 

The above results certainly deepen our theoretical understanding of four-fermion in- 
teraction models. They not only may give a definite theoretical constraint on structure 
of some phenomenological models p^] but also could provide a useful clue to the feature 
of these models at finite temperature and finite quark chemical potential, as was shown 
in the 3D GN model with A^^^ = 3. 



[1] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; 124 (1961) 246. 
[2] D.J. Gross and A. Neveu, Phys. Rev. D 10 (1974) 3235. 
[3] J. Goldstone, Nuovo Cimento 19 (1961) 154 . 

[4] D. Weingarten, Phys. Rev. Lett. 51 (1983) 1830; S. Nussinov, Phys. Rev. Lett. 51 (1983) 
2081; E. Witten, Phys. Rev. Lett. 51 (1983) 2351; C. Vafa and E. Witten, Nucl. Phys. 
B234 (1984) 173. 

[5] B. R. Zhou, Phys. Lett. B215 (1988) 364; B235 (1990) 412(E); Zhou Bang-Rong, Commun. 

Theor. Phys. (Beijing, China) 15 (1991) 319; 18 (1992) 347; 19 (1993) 377 . 
[6] V. A. Miransky, Dynamical symmetry breaking in Quantum Field Theory, World Scientific, 

Singapore (1993) and the references therein. 
[7] J. C. Colhns and M. J. Perry Phys. Rev. Lett. 34 (1975) 1353; M. Afford, K. Rajagopal, 

and F. Wiiczek, Phys. Lett. B422 (1998) 427; Nucf Phys. B537 (1999) 443; R. Rapp, T. 

Schafer, E. V. Shuryak, and M.Vefkovsky, Phys. Rev. Lett. 81 (1998) 53; J. Berges and K. 

Rajagopal, Nucl. Phys. B538 (1999) 215 . 
[8] M. Bubalfa, Phys. Rep. 407 (2005) 205. 



15 



[9] D. Diakonov, H. Forkel, and M. Lutz, Phys. Lett.B 373 (1996) 147. 

[10] B. Vanderheyden and A.D. Jackson, Phys. Rev. D61 (2000) 076004. 

[11] S. B. Rxister and D. H. Rischke, Phys. Rev. D69 (2004) 045011. 

[12] Zhou Bang-Rong, Commun. Theor. Phys. (Beijing, China) 47 (2007) 95. 

[13] Zhou Bang-Rong, Commun. Theor. Phys. (Beijing, China) 47 (2007) 520. 

[14] Zhou Bang-Rong, Commun. Theor. Phys. (Beijing, China) 47 (2007) 695. 

[15] H. Kohyama, Phys. Rev. D77 (2008) 045016. 

[16] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17 (1966) 1133; S. Coleman, Commun. 
Math. Phys. 31 (1973) 259. 



